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SUMMARY 

A theoretical model of the Raman rotational spectrum of a gas mixture was formu- 
lated and applied to an air mixture of oxygen and nitrogen. Theoretical spectra were 
calculated as a function of temperature and compared directly to spectra obtained ex- 
perimentally. The comparison, made at temperatures within the range 243 to 313 K, 
showed the model to be accurate. An error analysis was then performed yielding theo- 
retical and experimental temperature measurement accuracies. The theoretical error 
predicted for the temperature range used was ±20 K (~7 percent), which compared well 
with experimental data. The best experimental error was ±13 K (~5 percent). It was 
concluded that, within the limits of this experiment and the preceding accuracies, the 
temperature of a gas mixture can be measured by using Raman scattering. 


INTRODUCTION 

Laser systems have found many applications in the atmospheric sciences as tools 
for the remote measurement of meteorological variables. Some of the variables now 
being measured on an experimental basis include aerosol content, atmospheric density, 
atmospheric composition, and pollutant concentrations (refs. 1 to 3). Atmospheric tem- 
perature is a variable that at present is not being measured remotely by laser tech- 
niques. 

Cooney suggested a method that might be used to implement this measurement 
(ref. 4). Because the intensity distribution of Raman- scattered light is a function of 
temperature, it is theoretically possible to measure the air temperature by analyzing 
its Raman spectrum. By using an optical radar (lidar) system such as is described in 
reference 5, two or more portions of the spectrum of a remote volume of air can be 



measured simultaneously. The ratio of these measurements can then be related to the 
temperature of the volume of air. (Ratios of spectral intensities are used so that con- 
sideration need not be given to fluctuations in the index of refraction, water vapor, and 
particle concentrations, etc. , in the intervening atmosphere. ) The feasibility of using 
this method was demonstrated by Salzman, Masica, and Coney (ref. 6). In that study, 
it was shown that the temperature of a single gas system could be measured by analyzing 
Raman- scattered laser light. To relate the temperature of a volume of air to its Raman 
spectrum, it was necessary to expand that analysis so that application could be made to 

o 

a mixture of gases. 

This investigation was conducted to formulate a theoretical expression representing 
the Raman spectrum of air (and thus determine its temperature dependence) and to con- 
firm the results by comparison with experimental data. This report presents a summary 
of this theoretical and experimental analysis. The primary consideration here is with 
the nitrogen-oxygen combination in air. Because the ratio of the cross sections for ro- 
tational Raman scattering of nitrogen and oxygen is an integral part of this analysis and 
because the available (theoretical) value appears to be too inaccurate, specific attention 
was given to its measurement. Experimental data are also presented showing the varia- 
tion in the spectrum of air as a function of temperature for a temperature range of 243 to 
313 K and a pressure of 1 atmosphere. The data were obtained on a Raman spectro- 
photometer. 


ANALYSIS OF RAMAN SPECTRAL INTENSITY 

When light is incident upon a molecular system, a quantum mechanical interaction 
takes place between a small portion of the light energy and the rotational states of the 
system. Raman scattering of the incident light can result. Because of the quantized 
nature of the rotational energy levels, the Raman-scattered light is shifted to the dis- 
crete wavenumbers y n ±A y-, where y n is the wavenumber of the incident light and Ay. 

u J u th J 

is the wavenumber associated with the j quantized level, y Q - Ayj represents an 

absorption of energy by the molecule (Stokes interaction), and i>q + Ay^ represents an 
emission of energy (anti-Stokes interaction). The rotational Raman spectrum is a meas- 
ure of the intensity of the scattered light as a function of wavenumber. The molecular 
system of interest in this study is air; consequently, the Raman spectrum produced by 
a gas mixture must be analyzed. 

The spectral intensity I(y) of a mixture of independent gases is the sum of the 
spectral intensities of the component spectra corrected for the relative concentrations 
and Raman- scattering cross sections of the components: 


2 


^ixture^ ^omp 1^ + ^comp 2^ + K l I comp 3^ + • • • 


( 1 ) 



where K, Kj, etc. , are correcting constants normalized to component 1. (All symbols 
are defined in appendix A. ) Therefore, the theoretical Raman spectrum of air can be 
obtained by an independent analysis of the spectra produced by each constituent of air. 
Since nitrogen and oxygen make up approximately 99 percent by volume of the contents of 
air, the following analysis is limited to their contributions. Equation (1) thus becomes 
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with K related to the ratios of the concentrations and cross sections of oxygen and ni- 
trogen. 

A typical experimental Raman spectrum of a single gas system, in this case nitro- 
gen, is shown in figure 1. Each line profile in such a spectrum of a single gas can be 
written in terms of a line shape function S(y, Vj)/Aj normalized to a relative peak height 
Ij (Aj is a normalizing constant). It is shown in reference 6 that, for an analysis of a 
single-gas -system spectra, it is sufficient to consider only Ij (i. e. , one can approxi- 
mate the spectrum in fig. 1 by a line spectrum) and that, therefore, only peak heights 
need be measured in most experiments of this type. This approximation is not satisfac- 
tory for analyses of the oxygen-nitrogen system being considered herein because the line 
positions of the component spectra are different and partial or complete overlapping is 
likely. Since peak height is again generally the experimentally measured parameter, it 
is necessary to take into account contributions to each peak made by adjacent peaks. 

This can only be done by retaining in the theoretical analysis the line shape function 
S{v, i/j). The spectral intensity of an isolated peak is therefore given by the equation 




A i 


(3) 


where Aj normalizes the shape function to a peak value of 1, ^ is the wavenumber of 
the center of the peak, and j is the quantum number associated with the energy level 
involved in the Raman interaction producing the peak. The normalized shape function, 
Ij, and ^ are designated in figure 1 for the Stokes j -value of 12. The total spectrum 
of a single gas system is the sum of all the peaks resulting from each of the many pos- 
sible rotational energy level transitions: 


I(v) = 
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S(v, Vj) 
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(b) Anti-Stokes. 

Figure 1. - Experimental spectrum of a single gas system - nitrogen. (Mie and Rayliegh scattered light, i.e., Au-. - 0 is centered at 
20 492 cm' 1 .) J 



Substituting this equation into equation (2) yields the spectral intensity of air as a func- 
tion of wavenumber: 


‘air ' 1 ' 1 = 


V sN (* -i) 
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( 5 ) 


Raman Line Intensity 

An analysis of the relative peak intensities of an ideal line spectrum of a linear, 
diatomic molecule such as oxygen or nitrogen was performed in reference 6. Because 
these spectra result from molecular-energy-level transitions from the initial state j to 
j + 2 and j - 2, the line intensity Ij is determined by two component equations. It is 
shown in reference 6 that the line intensity of the Stokes half of the spectrum (j — j + 2) 
is 


where 


_ C(j + l)(j + 2) 9 
j 2j + 3 T 




A 


S 

Vi 



and of the anti -Stokes intensity (j - j - 2) is 
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where 
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(7) 


The constant C in equations (6) and (7) is a proportionality constant which includes a 
nuclear spin degeneracy factor, 9 is the rotational characteristic temperature of the 
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scattering gas, and T is the gas temperature. 

.It should be noted that 0 2 has a ground- state spin and therefore has a more com- 
plicated intensity distribution than is given by equations (6) and (7). However, because 
the variations due to this spin are only significant at very low temperatures («200 K), 
they were considered negligible for the purposes of this study. 

Equations (6) and (7) relate the spectral intensity of the ideal line spectrum of a 
single gas system to the gas temperature. The effect of this temperature dependence on 
the rotational envelope of the line spectrum of nitrogen is illustrated in figure 2. At the 
low temperature Tp most of the nitrogen molecules populate the low energy levels and 
the line spectrum is "grouped" near the ground state (j = 0). As the temperature of the 
gas is increased (e. g. , to Tg), more of the molecules populate the high energy states 
and the spectrum appears more evenly distributed. 

The alternating intensity of the odd and even j -lines in figure 2 is a result of the 
nuclear spin degeneracy of nitrogen. Because of the spin degeneracy, the intensity of 
the odd j -lines is one-half the intensity of the even j -lines. In the case of oxygen, the 
spin degeneracy is such that the even j -lines have zero intensity. 




Wavenumber 
(c) Temperature T3. 

Figure 2 . - Temperature variation of line spectrum ot single gas system - nitrogen. T j « T 2 « T 3. 
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Equations (6) and (7) are related to the relative line intensity factor 1^ used in equa- 
tion (5) by the equation 

for j — j + 2 (Stokes) 

for j — j - 2 (anti-Stokes) 

Rewriting equation (5) to reflect the distinction between portions of the Raman spectrum 
yields 
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t AS 
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Note that I? = I®(0 N ), etc. 

The same shape function is used for both sides of the spectrum since the same 
molecular levels are involved in the transitions producing the two sides. The difference 
is only in the direction of the transition; that is, the Stokes side results from an absorp- 
tion of energy by the molecule while the anti -Stokes side results from an emission of 
energy. 


Spectral Line Broadening 


An evaluation of the shape function involves an analysis of the effects of line broad- 
ening and instrument distortion on an ideal line spectrum such as that shown in figure 2. 
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Attention is first given to the profiles resulting from Doppler broadening, natural line 
broadening, and pressure broadening and to the profile resulting from their interaction. 
Instrument distortion is included in the analysis by applying the appropriate monochro- 
mator slit function to the resultant profile. Another factor classified as instrument dis- 
tortion is contributed by the detection and recording electronics. A theoretical analysis 
of this factor was not made. Rather it was assumed that bandwidth only would be affected 
and that the magnitude of this effect could be determined experimentally. 

It was expected that for this study the slit function would be the dominant factor in 
establishing the shape and bandwidth of the final theoretical spectrum. However, the 
other broadening factors were included in the analysis to determine the magnitude of 
their effect and to ensure completeness. 

Doppler broadening . - For this study the incident light energy is that produced by a 
continuous-wave, gas laser. Because the laser emission source is a gas of highly mo- 
bile ions, the output of the laser is Doppler broadened. The distribution resulting from 
this broadening is Gaussian in wavenumber (ref. 7) as given by 


S laser (l '> = ex P 


2.77/ x2 
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( 10 ) 


where is the half-width of the laser emission and uq is the center emission line. 
When applied to the calculation of the Raman line shape, i/q becomes uy 

Natural line broadening . - The energy levels of the target molecules represent a 
second source of line broadening. This broadening takes several forms, included among 
which are natural line broadening and those forms collectively referred to as pressure 
broadening (e. g. , collision, Stark, and impact broadening). The first, natural line 
broadening, results because of the uncertainty associated with the value of the emission 
wavenumber v of a molecule radiating from an upper level to a lower level (ref. 8). 

By evaluating the probability distribution of the energy of the upper and lower levels, a 
distribution as a function of wavenumber is obtained 


S nat (l/) 
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where t- and j- 9 are the lifetimes of the j** 1 and (j + 2)^ levels. The bandwidth of 
natural broadening calculated from the Lorentz-type distribution is 
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which corresponds to the uncertainty in v j as dictated by the Heisenberg uncertainty 
principle. 

Pressure broadening . - The second type of broadening, pressure broadening, men- 
tioned in the preceding section refers to those mechanisms whose effect increases with 
pressure. As with natural broadening, these generally produce Lorentz-type distribu- 
tions 


®press^ 


Rp/2 



( 12 ) 


where Rp is the bandwidth due to pressure broadening. The combined effect of natural 
and pressure broadening is simply 


S Raman*^ 


V 2 


<» - 
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where 





+ Rp 


and is the bandwidth due to natural and pressure broadening (ref. 8). 

The Raman- scattered light emitted by the molecules is the result of the interaction 
between the Doppler -broadened incident light energy and the natural- and pressure- 
broadened molecular energy levels. The resulting distribution about an emission wave- 
number i/j is given by the relation 

s Voigt (,/ ’ "j* = S laser (,/)S Raman (l '’ ^ (14) 


or 
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where v' is a dummy variable of integration. 

Equation (15) is commonly referred to as the Voigt function (ref. 8) and represents 
the broadened Raman spectral output for the j** 1 rotational level transition. Since a 
closed-form solution to this equation is not available, the emission line profile was ob- 
tained numerically. Figure 3 shows the theoretical profile resulting from the interaction 
of incident light having a line width of 0. 166 cm - * with a Raman energy level having a 
line width of 0. 075 cm - *. As is indicated in the figure, the bandwidth of the resulting 
emission profile is 0. 212 cm - *. 



Wavenumber, cm 


Figure 3. - Numerical solution to Voigt function using 0. 166 cm' 1 for laser line width and 0. 075 cm' 1 for 
Raman line width. 
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Instrument distortion . - The instrumentation used to measure spectral intensity as 

a function of wavenumber was a single-pass monochromator used in conjunction with 

photon-counting detection electronics. (This equipment is discussed in appendix B. ) 

Most instrument distortion of the spectral line shapes is produced by the monochromator 

as a result of "nonlinear transmission characteristics" of the slits (ref. 9). Because 

of the finite width of the entrance and exit slits and because of the imperfect resolution 

of the grating, a band of wavenumbers is transmitted. The distribution as a function of 

wavenumber of this transmitted light is commonly referred to as the slit function S ( j/) . 

s 

Because the slit function has a major effect on the shape and bandwidth of the meas- 
ured spectrum, its incorporation into the theoretical calculations is essential. An 
evaluation of the transmission characteristics of the monochromator used in this study 
suggested that either of the two commonly used slit functions (refs. 9 to 11) might be 
applicable. The two functions are the triangular slit function, having a half-width equal 
to the geometric slit width R g of the monochromator, and the Gaussian function 


S s (u) = exp 


2. 77 



(u - u s ) 2 


(16) 


where v c is the wavenumber setting on the monochromator. Both of these functions 
s 

were used analytically to determine which best reproduced an experimental Raman line 
profile. For the spectral system under study here, the Gaussian slit function was found 
to be the better approximation. Consequently, this function was used in all further 
analysis. 

The interaction of a slit function with a line-broadened emission profile is illustrated 
in figure 4. The Raman emission line profile Sy oi g t (^, ja) is centered at ia and the 
slits are centered at v a . The intensity of the light transmitted by the slits is propor- 
tional to the area (shaded in fig. 4) under the emission profile adjusted by the slit func- 
tion. Written in general, by letting y - u, this is 
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Figure 4. - Interaction between slit function and Voigt function. 


Substituting equations (15) and (16) into this equation yields 



Figure 5 shows the solution of equation (17) for R^ = 0. 166 cm" 1 , Rj^ = 0. 075 cm" 1 , 
and a slit width of 1. 84 cm" 1 . Again, the solution was obtained numerically. The band- 
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Figure 5. - Line profile resulting from combination of slit function and 
Voigt profile for slit width of 1.84 cm' 1 . Laser line width, 0. 166 cm' 1 ; 
Raman line width, 0.075 cm' 1 . 


width of the resultant profile is 1. 89 cm~*. 

Interference filter transmission . - In the case where interference filters are used 
in place of a monochromator, the shape function is equal to the Voigt function. There- 
fore, equation (5) becomes 



S Voigt^> V 



A. 

J 


^Voigt^j’ 


The total transmission of a filter is 
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*total S F (v)\ ir (v)d u 

The transmission of the filter is thus related to the temperature of the scattering gas 
through the relative peak intensities Z and l9 contained in I a - r . ' 

Calculation of a Theoretical Raman Spectrum 

Incorporating the results from the preceding analysis of line intensity and line shape 
into equations (8) and (9) yields a set of equations for the theoretical Raman rotational 
spectrum of air. These equations are presented in table I. The equations in this table 
are referred to by Roman numerals. 


TABLE I. - EQUATIONS FOR THEORETICAL RAMAN ROTATIONAL SPECTRUM OF AIR 
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Calculating a theoretical spectrum of air by using these equations involved the deter- 
mination and substitution of all pertinent constants, that is, the rotational characteristic 
temperatures 0^ and 9q and the appropriate bandwidths, and included a comparison of 
the component equations (II) to (VI) with experimental data. 

Rotational characteristic temperature . - The generally accepted values for the ro- 
tational characteristic temperatures of nitrogen and oxygen are 2. 862 K (ref. 12) and 
2. 082 K (ref. 13), respectively. The value given for 0^ was shown in reference 6 to 
reproduce accurately the line positions of an experimental spectrum of nitrogen. How- 
ever, a similar comparison between predicted and experimental line positions for oxygen 
made in the study reported herein showed considerable misalinement. (The experimen- 
tal procedure and apparatus used is discussed in appendix B. ) These data indicated that 
0 q equaled 2. 06 K. This value also agreed very well with data presented by Stoicheff 
(ref. 12). By using those data, a value of 2.065 K was calculated. This latter value 
was used in all further calculations of the air spectrum. 

Bandwidths . - The constants to be substituted into the shape function include the 
slit width R the bandwidth of the laser R^, and the bandwidth of the Raman energy 
level Rr,. The value substituted for R is simply the slit-width setting on the mono- 
chromator used to obtain the experimental spectrum. A value of 0. 166 cm' 1 for R L 
was supplied in the manufacturer's data for the laser. 

The bandwidth of the Raman energy level R^ was evaluated through a review of 
references 14 and 15. Reference 14 provided a discussion of the bandwidth as a function 
of the quantum number j, and reference 15 listed R^ as a function of the gas species. 
Experimental data presented in reference 14 showed a variation in R^ with j such that 
R^(j=7) = 0. 08±0. 01 and Rjj(j=12) = 0. 061±0. 01 cm’l To determine the significance of 
this variation on the Raman emission profile and on the spectrum profile, equations (15) 
and (17) were evaluated for R r of 0.085 cm 1 and 0.065 cm”*. Values used for R T 

_ 1 K _ i ^ 

and R were 0. 166 cm and 1. 84 cm , respectively. For this change in R r of 

S -j Xv 

0. 02 cm” , the bandwidth of Sy Qi t (y, ia) changed approximately 6 percent and the 

bandwidth of S(y, v-) changed approximately 1 percent. This variation was considered 
J _ i 

negligible for this study and it was assumed that R^ was a constant equal to 0. 075 cm . 

Because reference 15 indicated that the bandwidths of oxygen and nitrogen were the same, 

R^ was also set equal to 0. 075 cm" 1 . 

With the inclusion of these values for R g , R^, and R^ into equations (I) with 

K = 0, a theoretical line profile for a single gas system of nitrogen was calculated. A 

comparison was then possible between the theoretical and experimental bandwidths and 

line shapes. . The experimental bandwidth exceeded the theoretical by about 8 percent. 

This discrepancy was attributed to an additional distortion introduced by the response 

characteristics of the detection and recording electronics. The value of the slit width 

R used in the calculations was adjusted to account for this added distortion. This ad- 
s 

justment was made by first plotting the theoretical bandwidth as a function of R , as in 

s 
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3.2, 



Figure 6. - Line width of theoretical spectrum as function of slit-function half-width for 
laser line width of 0. 166 cm'l and Raman line width of 0.075 cm'^. 


figure 6. This plot was then used by measuring the bandwidth of the experimental spec- 
trum, locating this value on the abscissa, and using the corresponding value for R„ 

s 

(ordinate) in equations (I). 

With the inclusion of the broadening effects produced by the electronics, the theo- 
retical spectrum of nitrogen reproduces very well the experimental spectrum. Figure 7 
shows a comparison of a theoretical Raman spectrum of nitrogen with an experimental 
spectrum. The theoretical spectrum was calculated by using equations (I) with K = 0, 
and the constants and correction factor previously discussed. Peak heights were nor- 
malized according to the method described in appendix B. To proceed with the calcula- 
tion of a Raman spectrum of air then required only an evaluation of the constant K. 

Relative cross section . - As stated previously, the constant K is related to the 
Raman rotational cross sections o of nitrogen and oxygen and to the molecular concen- 
trations in air p of the two components. K is given by the equation 


K = ^0 

Vn 


(18) 


The percent concentrations of nitrogen and oxygen in air are well documented. Values 
obtained from reference 16 yield a molecular concentration ratio Pq/p n of 0. 27 at 
1 atmosphere. Apparently, no experimental data presently exist concerning the rota- 
tional cross sections of nitrogen and oxygen. The only analytical values found were 
those presented by Melfi (ref. 1). They are ct n = 3. 2x10 cm sr and 
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Figure 7. - Comparison of theoretical and experimental rotational spectra of nitrogen. 
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Oq = 5. 1x10 cm sr . Substituting these concentrations and cross sections into 
equation (19) resulted in a K-value of 0. 43. 

By using this value of K, a theoretical spectrum of air was calculated. A compar- 
ison was then made between the ratios of the oxygen to nitrogen peak heights measured 
from this spectrum and those measured from an experimental spectrum. Because these 
ratios did not agree, it was concluded that the accuracy of the theoretical value given 
for (1. 59) was inadequate. Consequently, data were obtained for the purpose of 

experimentally evaluating the relative cross sections of oxygen and nitrogen. 


RESULTS AND DISCUSSION 
Measurement of Relative Cross Section 

The relative cross section of oxygen and nitrogen was measured by comparing 
theoretical rotational Raman spectra of air to experimental spectra and determining that 
value of K which resulted in the best agreement. For this purpose, experimental spec- 
tra were obtained (by using the apparatus and procedures described in appendix B) for 
monochromator slit-width settings of 0.92, 1. 84, and 2.76 cm - *. Peak intensities of 
each of these spectra were measured and documented and equations (I) were appropri- 
ately normalized. 

Theoretical spectra were then calculated for various values of K at each slit- 
width setting. Since the relative molecular concentration of oxygen and nitrogen was 
assumed constant and fixed at 0. 27, this variation in K resulted directly in a variation 
of the relative cross section ct q/<? n - Theoretical spectra were calculated for ^q/^ n 
values of 1. 7, 2. 0, and 2. 3 at slit-width settings of 0. 92 and 1. 84 cm - * and 
values of 1. 7, 2. 15, and 2. 6 at a setting of 2. 76 cm - *. The theoretical intensity of each 
of the oxygen lines was then measured and plotted as a function of Oq/ct^. Figure 8 is 
typical of the resulting graphs. Shown are the variations obtained by using a slit width 
of 1. 84 cm - *. (A line designation of Og refers to the resolved j=3 oxygen line while 
a designation of Og-Ng refers to the line resulting from some degree of irresolution 
between the j=5 oxygen and j=3 nitrogen lines. ) 

By using these graphs the relative cross sections corresponding to each of the ex- 
perimental line intensities of oxygen were determined. An average was then calculated 
for each slit-width setting. Table II lists the lines used and the relative cross sections 
measured for the slit-width setting of 1. 84 cm - *. Average values obtained for settings 
of 0. 92, 1.84, and 2. 76 cm - * are 2. 11±0. 38, 2. 16±0. 35, and 2. 41±0. 27, respectively. 
The apparent discrepancy in these values is not attributable to a poor fit between theo- 
retical and experiment spectra, since in all cases they compared very well. However, 


18 



14 



Ol I I I I I I 
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Relative cross section, Oo/°N 

Figure 8. - Theoretical Stokes line intensity as function of relative cross 
section for slit width of 1. 84 cm' 1 . 


TABLE II. - OXYGEN EXPERIMENTAL LINE INTENSITY AND 
CORRESPONDING RELATIVE CROSS SECTION FOR 


SLIT WIDTH OF 1. 84 cm 



Stokes 




Anti-Stokes 


Line 

Intensity, 
counts/ sec 

Relative 

cross 

section 


Line 

Intensity, 

counts/sec 

Relative 

cross 

section 

°1 

33 

1. 45 


°3 

37 

1.70 

°3 

96 

2. 18 


°5 

85 

2. 17 

°5- N 3 

129 

2.24 


°7- N 5 

110 

2. 23 

°7 

121 

2.21 


°9 

95 

1. 85 

°9' N 6 

253 

2.02 


°ir N 8 

208 

2. 37 

°11 

106 

2.22 


°13 

83 

2. 38 

°13' N 9 

138 

1. 78 


°15- N 11 

104 

2.09 

°15 

61 

2.04 


°17 

55 

2. 75 

°19 

33 

2. 30 


°21 

23 

2.68 

°21 

18 

1.75 


°23 

15 

2. 81 
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the poor resolution resulting from a slit-width setting of 2. 76 cm~* limited the number 
of resolved oxygen lines and made base line and background intensities difficult to es- 
tablish. Consequently, confidence was placed in the results obtained from the better 
resolved data, that is, those obtained using slit widths of 0. 92 and 1. 84 cm - *. Taking 
the average of these two values resulted in a relative cross section of 2. 14±0. 36. Fig- 
ure 9 shows a comparison between the theoretical spectrum calculated using this value 
and an experimental spectrum. In this case the slit width was 1. 84 cm’l 


Measurement of Air Temperature 


The relations yielding the theoretical Raman rotational spectra of air are functions 

S AS 

of the air temperature through the line intensity factors Ij and Ij . The measurement 
of the gas temperature is accomplished by relating its temperature to the ratio of line 
intensities. Line ratios are used to eliminate extraneous effects such as encountered in 
the transmitting medium. This is especially important when considering a remote sens- 
ing system. 

Because of the complexity of equations (I), it is not possible to obtain directly a 
solution for temperature in terms of the line ratios. Consequently, the solution was ob- 
tained graphically. First, to reduce the number of possible line ratios to be considered 
to a manageable quantity, an error analysis was performed. Reference 6 showed that 
the standard deviation of the temperature was given by the equation 


AT 



where the slope is simply 


Slope = 


V T 1> V T 2> 
y T i> ’ yv 


T 1 " t 2 


(19) 


Therefore, theoretical spectra were calculated and all line intensities were measured for 
temperatures of 243 and 313 K. These are plotted in figure 10. To simplify notation, 
each peak in the spectra was assigned a number. The numbering started at the Rayleigh 
line and progressed arithmetically for each peak. For example, the Orj Stokes line 
became the S6 line, and the anti-Stokes line became AS9. This numbering system 
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(b) Anti-Stokes. 

Figure 10. - Comparison of theoretical spectrum of air for temperatures of 243 and 313 K. 









is illustrated in figure 10. By using equation (19) the standard deviation of the tempera- 
ture for theoretical peak intensities normalized to typical experimental values was cal- 
culated for all line-ratio combinations. This included Stokes-to-Stokes, anti-Stokes-to- 
anti-Stokes, and anti-Stokes-to-Stokes ratios. The best standard deviation predicted was 
±20 K. To determine the validity of equation (19), the standard deviation for each of 
these ratios was then compared to the standard deviations calculated from experimental 
data. To reduce the quantity of experimental data analyzed, attention was given only to 
those ratios which resulted in a theoretical AT of less than 25 K. Fifty-six line-ratio 
combinations satisfied this restriction. Because equation (19) is a conservative esti- 
mate of AT (ref. 6), it was expected that the experimental AT would be less than pre- 
dicted. 

Eight experimental spectra were obtained for comparison, two at each of the tem- 
peratures 243, 273, 299, and 313 K. The line intensities were then measured and ap- 
propriate ratio combinations calculated. Graphs were made showing the theoretical 
variation in these ratios as a function of temperature. Figure 11 presents a sampling 
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Figure 11. - Theoretical variation with temperature of ratios of Stokes 
lines 4 through 7 to Stokes lines 22 and 24. 
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of the resulting curves. The temperature corresponding to each experimental line ratio 
was determined by using these graphs, and experimental standard deviations were cal- 
culated. These deviations compared satisfactorily with those predicted from equa- 
tion (19). The ratios best used to measure the temperature of the scattering gas were 
chosen on the basis of the experimental standard deviations. Considering those ratios 
with the lowest deviations wherein any single line intensity is used only once resulted 
in the selection of the following ratios: S7/S24, AS5/AS18, and AS24/S4. The temper- 
atures obtained by using these ratios are compared to the actual air temperature in fig- 
ure 12. It is concluded that a measurement of the air temperature is feasible within the 
temperature range and predicted accuracy of ±20 K presented here. Through the use of 
the selected lines an experimental accuracy of ±13 K was obtained. Since the accuracy 
of this measurement is strongly dependent on statistics, it is assumed that increased 
spectral intensity, the use of the ratios of multiple lines, or the use of the ratios of 
spectral bands such as obtained with interference filters would improve that accuracy. 



Local gas temperature, K 

Figure 12. - Temperature measurement results calculated from the Stokes 
line ratio S7/S24, the anti-Stokes line ratio AS5/AS18, and the anti- 
Stokes/Stokes line ratio AS24/S4. 
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SUMMARY OF RESULTS 


The purpose of the study reported herein was to formulate a theoretical expression 
that represented the Raman rotational spectrum of air and to confirm the results by 
comparison with experimental data. Specific attention was given to the dependence of 
this spectrum on temperature. First, an expression representing the spectrum of a 
single gas system (nitrogen) was formulated and verified by comparing a theoretical 
spectrum to a rotational spectrum obtained experimentally. Then considering only the 
major components of air (nitrogen and oxygen), their concentrations, their molecular 
properties, and the ratio of their Raman cross sections, two single gas spectra were 
combined to produce one theoretical spectrum of air. Again, experimental data and 
theory were compared with good agreement. The following specific results were ob- 
tained: 

1. The measured Raman rotational spectrum of air results from the interaction of 
the Doppler-broadened laser output, the pressure-broadened Raman energy levels, and 
instrument distortion introduced by the monochromator and the detection and recording 
electronics. Table I presents the theoretical equations for the spectral intensity as a 
function of wavenumber. 

2. The ratio of the rotational Raman cross section of oxygen to nitrogen is 
2. 14±0. 36. 

3. The Raman rotational spectrum can be used to measure the temperature of air. 
The data obtained in this experiment yielded temperature measurement accuracies in 
agreement with those predicted for a single ratio by the statistics of the process, that 
is, approximately ±20 K. The use of selected experimental lines resulted in an experi- 
mental accuracy of ±13 K. By considering the average of several ratios and by increas- 
ing the statistics of the data, it is theoretically possible to improve this accuracy. Fur- 
ther improvement may be possible by considering the ratio of spectral bands such as 
are obtained by using interference filters rather than a monochromator. Because the 
basic theoretical spectrum is now available, it is a simple matter to superimpose upon 
it the transmission curves of a filter system and thus calculate the ratios of the trans- 
mission of those filters as a function of temperature as well as calculate temperature 
measurement accuracies. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 4, 1972, 

501-04. 
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APPENDIX A 


SYMBOLS 


A normalization constant 

C proportionality constant 

c speed of light, 3x10 ^ cm sec -1 

h Planck constant, 6. 626x10 J joule sec 
I intensity 

j rotational quantum number (initial state) 

K constant related to relative concentrations and cross sections of component gases 

_i 

k Boltzmann constant, 1. 380x10 joule K 
R half -width (FWHM), cm" 1 

5 shape function 

T temperature, K 

AT standard deviation of temperature, K 

6 rotational characteristic temperature, K 

v wavenumber, cm" 1 

vq wavenumber of incident light, cm” 1 
Ay wavenumber shift, cm” 1 

p molar concentration 

3 - 1 

a Raman-scattering cross section, cm sr 

r lifetime, sec 

Subscripts: 

F filter 

j quantum number of initial state 

L laser 

N nitrogen gas 

O oxygen gas 

P pressure 

R Raman 
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s 

Superscripts: 

anti- Stokes 

N nitrogen gas 

O oxygen gas 

c Stokes 
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APPENDIX B 


APPARATUS AND EXPERIMENTAL TECHNIQUES . 

Test Apparatus 

A laboratory Raman spectrophotometer was constructed to provide experimental 
data to compare with the theory. The instrument was used to make Raman spectrum 
measurements of gaseous dry air over a temperature range from 243 to 313 K at a pres- 
sure of 1 atmosphere. Basic units of the instrument are shown in the block diagram and 
accompanying photograph in figure 13. 

The light source used for excitation was a continuous -wave argon-ion laser which 
was tunable to a number of discrete wavelengths through the use of an intracavity Lit- 
trow prism. For this investigation, it was operated at a wavelength of 488 nm (20 492 
cm - *), which provided approximately 200 milliwatts of multimode power. The output 

Q 

beam was linearly polarized to better than 1 part per 10 with the electric vector 
oriented vertically. 

External optics (fig. 13) were employed to increase the scattered light radiance, as 
described in reference 17. After initially being directed toward the gas sample by a 
high -reflectivity mirror, the laser beam was passed through a crystalline half-wave 
plate which rotated the beam polarization by 90°. (This was done to minimize the 
Rayleigh- scattered light being scattered toward the collecting optics. ) The beam was 
condensed by a focusing lens to create a high-power -density scattering site in the gas 
sample. The optical path of the beam was terminated at a return mirror which refocused 
the beam back into the scattering site, thus creating repetitive passes of the beam 
through the gas sample and increasing scattered light radiance. The volume of maxi- 
mum light scattering or the "Raman light source" could be approximated by a horizontal 
cylinder 20 micrometers in diameter and 2 millimeters in length (ref. 17). Included in 
the external optics was a light sensor which monitored the power of a secondary laser 
beam. 

_3 

The gas samples were confined in a test cell approximately 10 cubic meter in 
volume. Temperature regulation of the gas was provided by flowing ethanol through the 
walls of the cell and using a controlled -temperature liquid bath. The temperature of 
both the cell and the gas was continuously monitored throughout each test. The optical 
windows in the cell through which the laser beam passed were parallel plates of coated 
quartz. At right angles to the axis of these windows were two optical ports which com- 
prised a portion of the collecting optics system. At gas temperatures below the ambient 
dewpoint, all windows and ports were continuously purged with dry nitrogen to prevent 
fogging. 
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Czerny-Turner grating 
monochromator 



Detection and recording electronics 



Figure 13. - Raman spectrophotometer. 
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Light scattered at right angles to the axis of the exciting laser beam was collected 
over a large solid angle by a lens which was a port in the cell. This lens, a field stop, 
and a secondary lens comprised an optical system (fig. 13) which magnified the image 
of the Raman source cylinder by a factor of 2. 5 and matched the optical acceptance angle 
of the monochromator. By rotating the magnified image through 90° with a dove prism, 
it was possible to efficiently fill the vertical entrance slit of the monochromator with 
collected light. The intensity of the collected light was enhanced by gathering that light 
which was scattered 180° away from the primary collection lens with a concave mirror. 

A crystalline quartz wedge depolarized the light before incidence on the entrance slit. 

Wavelength discrimination was provided by a single 3/4-meter Czerny-Turner grat- 
ing monochromator. Both slits were operated at widths of 0. 92, 1. 84, and 2. 76 cm'^ 
(mechanical slit widths of 20, 40, and 60 /urn) which provided both efficient coupling with 
the collecting optics and adequate wavelength resolution. For this investigation, the 
wavelength scanning speed was set at 0.2 nm per minute (~8.4 cm’%in). Because 
the monochromator was only a single-pass instrument, stray light was a problem if the 
intensity of the collected Rayleigh light and spurious reflected laser light was very high. 
Fortunately, the Rayleigh -scattered light from oxygen and nitrogen was highly polarized, 
and its collected intensity was significantly reduced by the 90° rotation of the exciting 
beam polarization. With this procedure, the collected light at the excitation wavelength 
was sufficiently reduced to keep stray light and primary grating ghosts at an insignificant 
level. 

Light emerging from the exit slit was focused down by a multielement lens onto a 
1 -millimeter -diameter, S-20 surface of a Bendix channeltron photomultiplier tube. A 
highly stable power supply provided the multiplier bias (normally 2250 V). The small 
size of the photocathode surface and the low-noise characteristics of the multiplier cir- 
cuit made it possible to operate the tube without cooling. When it was operated at room 
temperature, the combined noise and background count rate was about 20 counts per 
second. 

The signal from the photomultiplier tube was sent through a system of photon count- 
ing (or pulse counting) electronics for amplification and display. Pulse signals were 
first passed directly to a low-noise preamplifier for impedance matching and first-stage 
amplification before transmission to the remaining rack-mounted electronics. A pulse- 
shaping linear amplifier then conditioned the signal for input into a pulse-height analyzer. 
The pulse-height analyzer served as a discriminator-trigger circuit so that each signal 
pulse was weighted equally and converted into a narrow-width pulse. The output from 
the pulse-height analyzer was directed to a digital recording system. 

The digital recording system consisted of a 100-mHz scaler, a teletypewriter, and 
a paper -tape puncher. The scaler was gated to continuously count and record at 
1- second intervals. Its output was recorded on paper tape by the interfacing electronics 
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and the teletypewriter. The data on the paper tape were then read and plotted by a high- 
speed computer. The result was a digital Raman rotational spectrum such as is shown 
in figure 1. 


Experimental Procedure 

Once the optical components were alined and the detecting electronics were condi- 
tioned for optimum pulse counting, the operation of the equipment was essentially rou- 
tine, requiring only minor fine-tuning adjustments. Prior to each test, the electronic 
equipment and the laser were operated for a 30-minute warmup period. This was neces- 
sary for both stabilization of laser power and reduction of spurious counting noise. Dur- 
ing this period, the gas cell was purged and filled with the desired gas and its tempera- 
ture was adjusted. 

With the system filled and stabilized, the optical components were adjusted for an 
optimum Raman signal. The rotational Raman spectra were then obtained by scanning 
from the oxygen anti-Stokes j=23 line at 484. 7 nm (20 631 cm"^) upward in wavelength 
through the oxygen Stokes j=23 line at 491. 4 nm (20 350 cm”*). All tests were con- 
ducted in the dark to reduce the background count rate. During the tests, all system 
variables were continuously monitored and recorded. 

Because a direct comparison was to be made between the theoretical and experi- 
mental spectra, it was necessary to normalize the theoretical spectra to the correct 
maximum peak intensity. 

The normalization was accomplished by referring to equations (II). It is shown in 
reference 6 that these equations accurately approximate the peak intensities of the ni- 
trogen spectrum. Therefore, each experimental nitrogen line could be used to calculate 
a corresponding nitrogen j=6 peak intensity. Equations (I) could then be normalized 
such that the peak intensity of the j=6 line equaled the average of the intensities calcu- 
lated by using equations (II). 

To ensure accuracy, it was necessary to use only those lines that were either very 
nearly or completely resolved. These were determined by comparing nitrogen peak in- 
tensities obtained from a theoretical spectrum of air with those obtained from a line 
spectrum of nitrogen. Those that were the same or nearly the same were considered 
resolved. 

Before making the direct comparison, it was also necessary to scale the experimen- 
tal data for plotting purposes because the scan rate of the monochromator was not con- 
stant. Instead, it varied cyclically between 0. 00357 to 0. 00435 nm/sec (~0. 150 to 
~0. 182 cm” Vsec). Although small, this inconstancy caused a gross mispositioning of 
the spectral lines. To correct this difficulty, event markers were superimposed on the 
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spectra every 0. 5 nm (~21 cm”*) by the monochromator. Although of an insufficient 
frequency to correct for the variable scan rate, it did show that equations (VI) accura- 
tely located the nitrogen lines. The nitrogen lines, therefore, were used as scaling 
markers. By using equations (VI), the center of each even j-line was assigned its ap- 
propriate wavelength. The plotting routine used these positions to scale the experimen- 
tal data accurately. The normalization and scaling procedures described here enabled 
the direct and accurate comparison of theoretical and experimental spectra. 
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